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We study the phase diagram of strongly interacting matter including the inhomogeneous 
phase of chiral density waves (CDW) within the Polyakov loop extended Nambu - Jona- 
Lasinio (PNJL) model. We discuss the phase structure taking into account density and 
flavour dependence of the Polyakov loop potential parameter and temperature dependence 
of the four-point coupling constant of the NJL model. It is shown that the CDW phase ex- 
ists and that can be interpreted as a special realisation of quarkyonic matter. This fact is of 
particular interest because the existence of homogeneous quarkyonic matter is strongly con- 
strained. This also indicates that the study of inhomogeneous phases at finite temperatures 
and baryon densities are of special importance. 

I. INTRODUCTION 

An understanding of structure of the phase diagram of strongly interacting matter at finite baryon 
density is still an open problem despite of the many efforts that have been devoted to its study since the 
very beginning of the QCD era. The lack of experimental data and reliable tools of the direct QCD based 
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Q\ calculations are the main culprits of this situation. Nevertheless, the progress is still possible and our 
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experience must rely on the interplay between general arguments, models calculations and QCD lattice 
results accessible at the low chemical potential. A very good example of such approach is the PNJL 
model [1] which combines the chiral and the deconfinement order parameters linked to QCD at finite 
temperature through the parameter fits to the lattice data t^iM- 

The new insight into the high baryon density domain based on large N c expansion has been proposed 
by McLerran and Pisarski ||5|]. In the limit of a large number of colours N c — > °° the QCD phase diagram 
simplifies substantially. The low baryon density and temperature region is occupied by the confined, 
chirally broken phase where the matter pressure scales with N c as P ~ 0(1). The high temperature 
phase is deconfined and dominated by gluons with the P ~ N% scaling. Finally, at the low temperature 
and high baryon density there is a confined phase which scales as P ~ N c . The latter phase has been 
named quarkyonic. 

It is very natural to interpret the quarkyonic phase within the PNJL model as deconfined and chirally 
restored phase. This is based on the PNJL picture of the dense matter which consists of a degenerate 
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Fermi sea of quarks with the colorless particle-hole excitations at the Fermi surface. Then the P ~ 

n 

N c scaling at high density naturally arises. It was shown in paper [6J] that at N c = 3 the quarkyonic 
phase slightly precedes the chiral transition at finite density and low temperature. The existence of the 
quarkyonic phase within the PNJL model was also studied in Ref. [|7j]. 

Discussions of the QCD phase diagram frequently consider only homogeneous phases. However, at 
the high baryon density one also has to take into account the possibility of the crystal structures. The 

i n n 

very well known examples are large N c QCD [18], skyrmion crystals fl9|], LOFF phases 111011 . Overhauser 
effect 111 ill and chiral density waves I112ll . The competition between chiral density waves and uniform 



colour superconductors were also discussed II 1 311 . The chiral spirals has been already considered in the 
quarkyonic phase II 1411 . In this letter we would like to consider the chiral density waves in quarkyonic 
matter from the point of view of the PNJL model. It is an interesting problem to check how the chiral 
density waves are influenced by the presence of the Polyakov loop field and what kind of the feedback 
it generates. We also pay a particular attention to the relation between the inhomogeneous chiral phase 
and the quarkyonic matter. 



II. CHIRAL DENSITY WAVES IN THE PNJL MODEL 

We consider the PNJL model with two light quarks and three colours with the parametrisation of the 
Polyakov loop based on paper [3] and the NJL part based on paper fl 1 311 



L = W^¥+ G[(W) 2 + (WsW) 2 ] - U(®, T) 



(1) 



where \\r is a massless quark field and D jJ = d jJ — iA^ is covariant derivative. The SU (3) gauge field 
= (Ao,0), Ao = gA C Q% a /2 and X a are the Gell-Mann matrices. The effective potential describing the 
traced Polyakov loop 
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in the fundamental representation takes the form y[] 



a(T) 



^a(7/)4> 2 + b(T) ln[l - 64> 2 + 8<£ 3 - 34> 4 ] 
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1.75 



(2) 



where A4 = iAq is a colour gauge field and the temperature Tq = 270 MeV describes the deconfinement 
transition in a pure gauge sector. In equation ([2b we already used the constraint that the expectation 



value of the Polyakov loops (<&), (<&*) are real 



160 which gives $ = <I>* at the mean field level |3|, m 



Above treatment may be improved with fluctuations that lead to a difference between the expectation 
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values of the traced Polyakov loops at non-zero baryon densities II 1711 . However, in the present work, as 
a practical procedure, we consistently keep 4> = Among the literature, one encounter also another 
approach. Fields and 4>* are treated as independent in the minimalization of the grand thermodynamic 
potential II 1811 . 

The coupling constant G = 5.024 GeV 2 and the three dimensional momentum cut-off A = 0.653 
GeV was used for regularization of the divergent vacuum contribution which reproduce the correct 
values of the pion decay constant and chiral condensate. 



We are working at the mean field level with the chiral density wave ansatz Ill2ll 



(\\n\f) = Mcosq-x, (\|/i'Y5'W) = Ms'mq-x, 



(3) 
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where the third direction in isospin space was chosen arbitrarily. The colour field in a Polyakov gauge 
creates a constant background field A4 = "k^ty which is related to the Polyakov loop through the equation 
4> = (l+2cos4>)/3. 

The standard procedure of the mean-field calculation in the Matsubara formalism leads to the ther- 
modynamic potential of the form 

" = ^ + ^*^)-il/(W Trtorl( *' p1 

where G)„ = (2n + 1)%T. The inverse propagator in Nambu-Gorkov space reads 

«Yo(oJn +M + i/j) +y(p - 5Y5T34) +M 

/Yo(<»h -A 4 - i/u) +y(p- jj 5 ^q) +M 

where /j is a quark chemical potential. After diagonalization of the propagator one finally arrive at the 
formulae 

& ^ M 2 M 2 F 2 q 2 ,„ f A d 3 p „ 

V y 1 ' AG 2M 2 J (2%y 

r-A ^3 



(5) 
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where Mq = 0.301 GeV is a constituent quark mass at zero temperature and density. The regularization 
through the 3-dim momentum cut-off was introduced after paper J13I1 . The energy eigenvalues are given 
by the expressions 



E± = \lp 2 +M 2 + ^-±\/(q-p) 2 +M 2 q 2 , E = ^p 2 +M 2 . 



Let us notice that potential © reduces to NJL model prediction in the deconfined limit 4> = 1 
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III. THE PHASE DIAGRAM 



The global minima of the potential © as a function of temperature and baryon density describe the 
phase diagram of the strongly interacting matter. This prescription leads to the self-consistent equations 

da _ da _ da 

dM 34> d\q\ 
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FIG. 1 : The phase diagram of the PNJL model with inhomogeneous chiral wave. The dashed curves are lines of 
the deconfinement phase transition. The left panel shows the diagram for the constant Polyakov loop potential 
parameter To = 270 MeV, whereas the right panel for the density dependent parameter Tq(jj) given by equation 
0. The gray colour describes the inhomogeneous CDW phase. The dotted curves are lines of the transition to 
quarkyonic phase. These lines essentially coincide with the lines of the first order phase transition to CDW phase 
at higher density. 

Fig. 1 shows the the PNJL model phase diagram including the inhomogeneous chiral density wave 
phase (CDW) marked with the gray colour. This phase is surrounded by the line of the first order 
phase transition where there are jumps in the values of the all order parameters M, and \q\. The 
jump in <I> is induced by the first order phase transition in M and it is not a mark of the transition to a 
deconfinement phase. One expects rather a crossover here and the exact place of this transition is not 
fixed unambiguously. We define the transition line as a place where the derivative dQ/dT reaches its 
maximum unless stated otherwise (see subsection MB). Using our definition the lines of deconfinement 
phase transition are shown in Fig. 1 as dashed curves. For the Polyakov loop potential parameter Tq = 
270 MeV the chiral and deconfinement transitions almost perfectly coincide at zero density, however, 
they split around the point (T,/j) = (0.21,0.12) GeV (left panel of Fig. 1). It is interesting to note that 
this is also the place where the quarkyonic phase appears. The boarder line for the quarkyonic matter is 
defined as a place where the value of the quark chemical potential exceeds the value of the constituent 
quark mass flg]. This line perfectly coincides with the first order phase transition line of the CDW phase 
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for /j > 0. 12 GeV. It happens because at the first order phase transition there is a large jump in the value 
of the constituent quark mass M which drops below the value of /j. One can conclude that within the 
PNJL model the quarkyonic matter can be treated as confined, however, spatially inhomogeneous phase. 

A. Deconfinement and the baryon density 

It is important to remember that the critical temperature of the deconfinement phase transition de- 
creases with increa sing baryon number density. This results in Tq parameter dependence on the quark 



chemical potential 111911 . It can be intuitively understand in the picture of overlaping hadrons. At zero 
density the finite temperature causes fluctuations of mesons and baryon-antibaryon pairs up to some 
critical value of T c . Above this temperature hadrons overlap with each other and there is a possibility 
of a colour flow in space which one can interpret as the process of deconfinement. At higher baryon 
chemical potential there are already some number of hadrons present in a medium. Then a new critical 
temperature, lower then T c , is sufficient to create appropriate number of hadrons which start to overlap. 
In paper I120ll K. Fukushima tried to estimate such dependence from the statistical model data which we 
also adopted here 

T (v)=T -9bv 2 , (7) 

where b = 1.39 • 10~ 4 MeV 1 and for Tq(/j) < 1 MeV the phase is defined as deconfined. On the right 
panel of Fig. 1 the dashed curve describes the deconfinement phase transition for the density dependent 
parameter O. As expected the line bends toward the density axis, but still, there is a place for the 
inhomogeneous phase in the quarkyonic region for temperatures below 0. 16 GeV and chemical potential 
above 0.22 GeV. This needs to be confronted with the statement that the existence of homogeneous 



quarkyonic phase is inconsistent with the prediction of statistical model |20|]. The introduction of non 



uniform phases open up a window for realisation of the quarkyonic matter (see also Fig. 4). 

B. Deconfinement and the number of flavours 

A value of the temperature parameter Tq depends also on the number of active flavours. For the model 
which contains two degenerate flavours one should reduce the temperature parameter to Tq = 208 MeV 
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2LI19H. The phase diagram for this parameter is shown on the left panel of Fig. 2. First of all one looses 
a good coincidence between the chiral and deconfinement phase transitions at zero density, nevertheless, 
the CDW phase remains intact. The confined (quarkyonic), inhomogeneous phase appears at chemical 
potential larger then /j = 0.21 GeV and temperatures below T = 0.16 GeV. For lower baryon densities 
and higher temperatures the CDW phase still exists but it is deconfined. 
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FIG. 2: The phase diagram of the PNJL model with inhomogeneous chiral wave for the Polyakov loop potential 
parameter To = 208 MeV which corresponds to Nf = 2 active degenerate flavours. On the left panel the NJL 
parameter is set to G = 5.024 GeV -2 . On the right panel the temperature dependence of the NJL parameter 
G(T) is taken into account according to equation ©. Behavior of the Polyakov loop is disturbed by the first 
order transitions into the inhomogeneous phase and consequently deconfmement transition is determined with 
uncertainty within the blue band. 

The phase diagram at To = 208 MeV is not satisfactory in a sense that at zero density there is a mis- 
match between the chiral and deconfinement phase transitions which contradicts the lattice results 112 ill . 
One can consider a possibility that this problem is a consequence of the wrong value of the coupling 
constant G at non-zero temperature. The effective coupling constant G is in principle a function of tem- 
perature and density. It contains contributions which follow from the integration of the gluonic degrees 
of freedom. At higher values of T gluons interaction with quarks weakens which in turn influences the 
effective four-quark interaction in the same way. Thus the coupling G decreases with increasing tem- 
perature which lowers the critical temperature of the chiral phase transition. If one attributes the whole 
mismatch between the chiral and deconfinement transitions at zero density to the wrong value of G then 
one can try, at least at the phenomenological level, to reestablish the agreement changing its value into 
the new one G(T C ) = G c in such a way, that both transitions has the same critical temperature T c . We 
assume here that the coupling G is a linear function of temperature 



G{T) = G(l - (T/T C )) + G C (T/T C ), 
G = 5.024 GeV 2 , G c = 4.221GeV 2 , 



(8) 



167 MeV. 



According to our fit the change in the constant G between zero temperature and T c is of the order of 15 
per cent which is not much. Then one can treat equation ([8]) as a series expansion in the temperature 
around T = 0. The density dependence of G is neglected since it is less important at large N c limit where 
quarks decoupled from the gluonic degrees of freedom. 
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The right panel of Fig. 2 describes the phase diagram where the NJL coupling constant G is given by 
the linear function from eq. ©. It is clearly seen that the chiral and deconfinement phase transitions 
coincides at zero density to a good approximation. Let us mention that the usual definition of the 
deconfinement transition line as a place where the derivative d<&/dT reaches its maximum is not quite 
useful in a situation where is a discontinues function of temperature and the points of discontinuity 
are close to the expected maximum value of d<&/dT . In such situations we define the deconfinement 
transition as a place where the Polyakov loop 4> takes the value 0.31 ±0.07. This is rather a modest 
value, however, such range is suggested when one considers the deconfinement phase transition for the 
homogeneous phases where the transition line is defined in a standard way. Using our definition the 
line of deconfinement phase transition is shown in Fig. 2 (later in Fig. 4) as a dark (blue) band. At 
low temperature and high baryon density one recovers the results of the NJL model with a constant G 
coupling. 
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FIG. 3: The chiral phases of the PNJL model with inhomogeneous chiral wave. The NJL model parameters are 
G = 5.024 GeV 2 and A = 0.653 GeV. The temperature dependent part of the potential © is also regularized by 
the finite cut-off. In such a regularization the CDW phase is limited to the high density region of the diagram. 



C. Triple point 



The concept of a triple point in the QCD phase diagram refers to various situations. For instance, 
the existence of a triple point between hadronic, color- superconducting and quark-gluon sectors was 
discussed in Ref. [22]. In the context of the present work, we refer to a very recent idea, namely, to the 
triple point in which hadronic, quarkyonic and quark-gluon phases meets together. Such a possibility was 
emphasized in paper ||23ll . Although in Figs. 1,2 the CDW phase exists even at zero density region (there 
is no space for triple point), but the location of the inhomogeneous phase depends on the regularization 
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parameter A. In Figs. 1,2 the cut-off parameter in the temperature dependent part of the potential © 
was sent to infinity. If one keeps this parameter at the constant value A = 0.653 GeV then the phase 
diagram changes and it is shown in Fig. 3. In this case the triple point actually appears at the point 
where the inhomogeneous phase ends (73,^3) = (0.22,0.2) GeV, what is in agreement with the results 



of Ref. 112411 . However, one has to remember that the location of this point is strongly dependent on the 
model parameters and the regularization method. 

For the quark chemical potential below /J3 the phase line describes the continuous phase transition 
whereas above /J3 the line is of the first order. When a non-zero current quark mass is turned on then the 
continuous phase transition changes into a smooth crossover and only the island of CDW phase remains 
on the diagram. 

IV. SUMMARY AND CONCLUSIONS 

We have discussed the phase diagram of the strongly interacting matter including the spatially de- 
pendent chiral density waves in the PNJL model. It was shown that the inhomogeneous CDW phase 
exists and dominates over the chirally restored phase in a large domain of the phase diagram. We have 
also pointed out that the chiral density wave can be interpreted as a special realisation of the quarkyonic 
matter. This is an interesting possibility, particularly, when the homogeneous quarkyonic phase was 



strongly constrained by the results of the statistical model H20I1 . Indeed in a case of two flavours the 
phase diagram with only homogeneous phases is shown in Fig. 4 (left panel) where the effects of flavour 
and density dependence of the Polyakov loop potential parameter To(Nf,fi) as well as the temperature 
dependence of the NJL coupling constant G(T) were taken into account. This last dependence reflects 
the fact that at higher temperature the four-point quark interaction should weaken which in turn let the 
chiral and the deconfinement transitions stay at the same critical temperature at zero baryon density. 
In the diagram the deconfinement phase transition precedes the chiral phase transition almost in all the 
domain. Only at high density and low temperature the transitions start to coincide. In such situation 
there is no space left for quarkyonic matter in accordance with the conclusion of paper ll2oll . 

The phase diagram in the right panel of Fig. 4 sums up all physical effects that we discussed in the 
previous sections and is the best candidate for the QCD phase diagram. The left panel of Fig. 4 should 
be compared with the right panel of the same figure. Let us notice a clear change in the order of the phase 
transitions. The phase transition to the CDW phase precedes the deconfinement phase transition. This 
fact also opens a window for the quarkyonic matter. The transition into the quarkyonic phase is close 
to or coincides with the first order phase transition into the CDW phase. Then the matter is spatially 
inhomogeneous but still confined. Let us remind that the details of the phase diagram depends on the 
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FIG. 4: Left panel: the phase diagram of the PNJL model with homogeneous phases only. Above the point (T, fi) 
= (0.085, 0.263), the nature of the chiral transition changes from first to second oreder. Right panel: the phase 
diagram of the PNJL model with inhomogeneous chiral density wave. The Polyakov loop potential parameter 
Tq(p) is given by © with Tq = 208 MeV. The finite 3-dim momentum cut-off A = 0.653 GeV regularized the 
potential © and the NJL coupling constant G(T) depends on temperature through the eq. ([8]). The grey colour 
depicts the CDW phase. The dark (blue) band shows a location of the deconfinement phase transition. 

temperature parametrisation of the NJL coupling constant G(T). Nevertheless, at low temperatures one 
should expect the same pattern of phase transitions as given in Fig. 4. 

For the model parameters we choose the CDW phase exists up to a zero density line. However, the 
low density and high temperature region is strongly affected by the temperature fluctuations which are 
neglected in the mean field approximation. These fluctuations probably melt the "crystal" structure of 
the CDW phase and the triple point is expected to appear on the phase diagram eventually. 

It is an interesting task for the future work to compare the phases of the chiral density waves and the 



chiral spirals [14]. However, it requires the implementation of the chiral spirals within the NJL model 
in the first place. This is a subject of certain importance because the quarkyonic matter, if exists, would 
be most probably of an inhomogeneous nature. One should also check the dependence of the results 



against the regularization scheme and the influence of the finite current quark mass II25I1 . Finally, it is of 
great interest to consider more general ansatz then © to study the possibility of the creation of different 
crystal structures. 
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